CRITICAL FLAVOR NUMBER IN THE THREE DIMENSIONAL THIRRING 

MODEL 
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We present results of a Monte Carlo simulation of the three dimensional Thirring model with the 
number of fermion flavors Nf varied between 2 and 18. By identifying the lattice coupling at which 
the chiral condensate peaks, simulations are be performed at couplings g 2 {Nf) corresponding to 
the strong coupling limit of the continuum theory. The chiral symmetry restoring phase transition 
is studied as Nf is increased, and the critical number of flavors estimated as Nf c — 6.6(1). The 
critical exponents measured at the transition do not agree with self-consistent solutions of the 
Schwinger-Dyson equations; in particular there is no evidence for the transition being of infinite 
order. Implications for the critical flavor number in QED3 are briefly discussed. 
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The study of quantum field theories in which the 
ground state shows a sensitivity to the number of fermion 
flavors Nf is intrinsically interesting. According to cer- 
tain approximate solutions of Schwinger-Dyson equations 
(SDEs), in d = 3 spacctime dimensions both quantum 
electrodynamics (QED3) and the Thirring model display 
this phenomenon. Both models have been proposed as ef- 
fective theories describing different regions of the cuprate 
phase diagram, Thirring describing the superconduct- 
ing phase, while QED3 supposedly describes the non- 
superconducting ''pseudogap" behaviour seen in the un- 
derdoped regime [l|, [2| • The Thirring model is a theory 
of fermions interacting via a current contact interaction: 
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where ipi,ipi are four-component spinors, to is a parity- 
conserving bare mass, and the index i runs over Nf dis- 
tinct fermion flavors. In the chiral limit to — > the La- 
grangian (JTJ shares the same global U(l) chiral symmetry 
if) 1 ► e M75 -0, ip 1— > -0e M75 as QED3. Since the coupling 
g 2 has mass dimension — 1, naive power-counting suggests 
that the model is non-renormalisable. However [1,H, @, 
an expansion in powers of l/Nf, rather than g 2 , is ex- 
actly renormalisable and the model has a well-defined 
continuum limit corresponding to a UV-stable fixed point 
of the renormalization group (RG). The l/Nf expansion 
may not, however, describe the true behaviour of the 
model at small Nf. Chiral symmetry breaking, signalled 
by a condensate (iptp) ^ 0, is forbidden at all orders in 
l/Nf, and yet is predicted by a self-consistent SDE a 
proach 0, [1, 0, Hi ■ SDEs solved in the limit g 2 -> 00 
show that chiral symmetry is spontaneously broken for 
Nf < Nf c ~ 4.32, close to certain predictions of Nf c for 
non-trivial IR behaviour in QED3 Q. Based on these 
results, at Nf — Nf c the model is expected to undergo 
an infinite order or conformal phase transition, originally 
discussed by Miranskii et al in the context of quenched 



QED4 [Ut| ■ Using different sequences of truncation, 
however, other SDE approaches have found values of Nf c 
ranging between 3.24 [j] and 00 [||. 

Since for Nf < Nf c the chiral symmetry breaking tran- 
sition corresponds to a UV-stable RG fixed point at a 
critical g 2 , which can be analysed using finite volume 
scaling techniques, study using lattice Monte Carlo sim- 
ulation has proved possible with relatively modest com- 
puter resources [T2I Il3l . [lH Il5j . It has been shown that 
for Nf < 5 the model is in the chirally broken phase at 
strong enough coupling g 2 < 00, providing a lower bound 
on Nfc, and the critical exponents found to vary with 
Nf. This should be contrasted with numerical studies 
of QED3: in this case the phase transition in the con- 
tinuum limit at Nf = Nf c is an IR-stable fixed point 
(see, eg. [l|| H3): and SDE predictions for the separa- 
tion of scales parametrised by the dimensionless quantity 
{ijjip} /a 2 are so small that enormous lattices are required 
to establish whether the model with a particular Nf lies 
in the broken or symmetric phase (l8j . A determination 
of Nfc for QED3 (whose value may have profound im- 
plications for the cuprate phase diagram 2]) by purely 
numerical means currently appears very difficult. 

As already noted, Thirring and QED3 share the same 
global symmetries; moreover in the strong coupling limit 
the l/Nf expansion predicts the existence of a massless 
spin-1 i/jip bound state in the Thirring spectrum [5], 0]. 
The RG fixed point of both models is therefore a theory 
of massless fermions interacting via massless vector ex- 
change between conserved currents. It is therefore plau- 
sible that the two models have the same Nf c , and that 
simulation of the Thirring model is the most effective 
means to determine its value. Here, we present results 
based on numerical simulations with Nf = 2, . . . , 18 in 
an effort to determine Nf c with unprecedented precision, 
exploiting a strategy of identifying the value of the lattice 
coupling g 2 corresponding to the strong coupling limit of 
the continuum theory. 

The lattice action we have used is based on [l§] and is 



2 



as follows: 
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where x, x are staggered fermion fields and the flavor 
index i runs over TV species. We have introduced an aux- 
iliary real- valued link field Ay, superficially resembling a 
photon field, so that the lattice action can be expressed 
as a bilinear form. The formulation is not unique; Gaus- 
sian integration over the Ay in ([2]) results in a lattice 
action resembling the continuum form ([1]) in that all in- 
teractions remain of the form XXXX f° r arbitrary N, but 
some of which when re-expressed in terms of continuum- 
like Dirac spinors are of non-covariant form. As argued 
in [13j ]. in the 1/Nf expansion these unwanted contri- 
butions are probably irrelevant, but more care may be 
needed when discussing the UV fixed-point theory. 

For d — 3 N staggered fermion species describe Nf — 
2N continuum flavors. We employed the Hybrid Monte 
Carlo (HMC) algorithm to simulate even Nf [HI and the 
Hybrid Molecular Dynamics (HMD) algorithm for odd or 
non-integer Nf [l~5| . For the HMD simulations we used a 
small enough fictitious time step At < 0.0025 to ensure 
that the 0(N 2 At 2 ) systematic errors in the molecular 
dynamics steps were smaller than the statistical errors of 
the various observables. 





FIG. 1: Leading order vacuum polarisation in lattice QED. 

Next we review the discussion of [T^ regarding the 
non-conservation of the vector current in the lattice 
Thirring model. The leading order 1/N quantum cor- 
rections to the photon propagator in lattice QED3 are 
sketched in Fig. [TJ The second diagram, arising from the 
gauge- invariant form X x e ""Xx+fn is peculiar to the lat- 
tice regularisation, but is required to ensure transversity 
of the vacuum polarisation tensor, ie: 



E FWa:) - IWi - /0] =0. 



(3) 



For the Ay propagator in the lattice Thirring model @ 
however, the second diagram is absent, and the transver- 
sity condition ([3]) violated by a term of 0(a _1 ). Transver- 
sity of n M „ is crucial to the renormalisability of the 1/Nf 
expansion; fortunately, the impact of the extra diver- 
gence can be absorbed by a wavefunction renormalisation 
of Ay, and a coupling constant renormalisation 
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where J(m) is the value of the integral contributed by 
the second diagram in Fig. [TJ The physics described by 
continuum 1 /Nf perturbation theory occurs for the range 
of couplings g 2 R £ [0, 00), ie. for g 2 € [0, gj- m ) where to 
leading order in 1/N g 2 im = § for m = 0. The strong 
coupling limit is therefore recovered at g 2 — g 2 m . For 
g 2 > g 2 m the auxiliary propagator becomes negative, and 
the lattice model no longer describes a unitary theory. 

As already discussed, chiral symmetry breaking is ab- 
sent in large-TV/ calculations. It may well be that the 
value of the second diagram in Fig. [TJ is considerably al- 
tered in the chirally broken vacuum expected for Nf < 
Nf c . In this study we use lattice simulations to estimate 
the value of g 2 m as a function of Nf . Fig. [2] shows (^ip) 
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FIG. 2: Chiral condensate (ipip) vs. g~ 2 . N f = 6, m = 0.01 
on a 16 3 lattice unless otherwise stated. The dashed line 
denotes a fit of the 16 3 , Nf = 6, m = 0.01 data to an EoS of 
the form ([5]) taking finite volume scaling into account. 

data as a function of g~ 2 for various Nf, m and lat- 
tice volumes. The condensate is non-monotonic, show- 
ing a clear maximum for Nf = 6 flavors at g~ 2 ~ 0.3. 
The peak position at g 2 = j^, unlike the value of the 
condensate itself, is independent of both lattice volume 
and fermion mass, indicating that its origin is at the UV 
scale. Since in an orthodox description of chiral sym- 
metry breaking |(?/>'0)l is expected to increase with the 
strength of the interaction, we interpret the peak as the 
point where unitarity violation sets in, ie. g^ax ~ Sum' 
and hence identify the peak with the location of the 
strong coupling limit. 

The figure also shows that g^ax depends on Nf (Cf. 
Fig. 3 of Ref. [3). For Nf < N fc , it is possible to fit 
condensate data taken at q 2 < .g^ ax to an equation of 
state (EoS) of the form p]f 



m = A(g- 2 - g- 2 )(iji') p + B(^) 6 



(5) 



which describes a continuous transition in the limit m — > 
at g 2 — g 2 to a chirally symmetric phase, characterised 
by critical exponents S and a "magnetic" exponent = 
(S—p)^ 1 . Ona dataset taken onl2 3 ,...,32 3 taking finite 
volume scaling into account we have found the best fit 
for Nf — 6 given by g~ 2 = 0.316(1), S = 5.75(13), p = 
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1.18(2), to be compared with corresponding quantities 
tabulated for smaller Nf in [la ]. Since g 2 $ 3^ ax , the 
EoS fit for Nf = 6, shown for the 16 3 data by a dashed 
line in Fig. [2J is of borderline credibility. 

Fig. [2] also shows that for Nf — 8 the peak moves to 
the right, and is considerably less pronounced. In Fig. [3] 
we plot g-l, vs. N f for N f G [2, 18]. For N f < 6, g^L 
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decreases smoothly and monotonically, but in the interval 
between 6.5 and 6.6 there is a sudden sharp increase; for 
6.6 < Nf < 12, g~ 2 ~ 0.43(2) is roughly constant, and 
then increases for Nf still larger. With the identification 
of g 2 layL with the strong coupling limit of the continuum 
theory, as argued above, it is natural to ask whether the 
transition at Nf ~ 6.5 has any correlation with the chiral 
order parameter. Fig. [4] shows (V'V'(ffmax)} vs - Nf. This 
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FIG. 4: vs. Nf for m = 0.01. Lines denote EoS 

fits discussed in the text. 

plot necessarily puts an upper bound on the condensate. 
It is difficult to draw any conclusion other than there 
being a chiral symmetry restoring phase transition at 



N fc = 6.6(1). 



(6) 



shown in Fig. [3] and the predictions of the 1/Nt expan- 
sion in the regime Nf > Nf c . In particular, there re- 
mains a disparity between g^^Nf) and the large- Nf 
prediction g 2 m = |, which may be due to subleading 
corrections in l/Nf, or to large finite- volume corrections 
described by the conformal field theory expected in the 
limit Nf -> Nf c+ . 

Fig- His interesting because for the first time it presents 
lattice data in a form suitable for direct comparison 
with SDE predictions. One such study of chiral sym- 
metry breaking in the d-dimcnsional Thirring model 
with d 6 (2, 4) using this approach was by Itoh et al 
0- They calculated the dressed fermion propagator 
S(p) = [A(p 2 )]/) + B(p 2 )}^ 1 , by first exploiting a hid- 
den local symmetry to fix a gauge in which A = 1, 
and then approximating the auxiliary propagator and 
fermion-auxiliary vertex by their forms to leading or- 
der in l/Nf. This enabled a solution for the self-energy 
function = B(x)/A, with A a UV cutoff and 

x = (p/A) d ~ 2 , in the strong coupling limit g 2 — 
yielding a dynamically-generated fermion mass M: 
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where N fc (d = 3) = 128/3tt 2 ~ 4.32. The chiral order 
parameter follows via the relation {tptp) oc S'(a; = 1), 
leading to the strong coupling prediction 



(tpip) oc AtM 1 . 



(8) 



Spontaneous chiral symmetry breaking occurs for Nt < 
Nf c in qualitative, but not quantitative agreement 
with Fig. 2J Eqn. predicts an infinite-order phase 
transition. Its nature is further elucidated using the 
anomalous scaling dimension of the ipip bilinear = 
d]xi{ipi/})/d]iiA and the relations for the critical expo- 
nents 77 and 5 [l9| : 
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From ([5]) we obtain 7^^ = (d — 2)/2 and hence 77 = 2, 
5 = 1. This scenario has been termed a "conformal phase 
transition" [ll[, and has also been exposed by SDE ap- 
proaches in QED3 [l6| and quenched QED4 [lO] , the con- 
trol parameter being respectively Nf and the fine struc- 
ture constant a. 

The SDE analysis of was later extended to cover 
g 2 < 00 by Sugiura [|[. In this case the chiral transition 
takes place for Nj < Nf c , and the solution for the order 
parameter is of the form 



(ipip) oc A d - 2 M. 



(10) 



Eqn. © is the main result of this Letter. 

In order to refine this picture, further theoretical anal- 
ysis is needed to establish contact between the results 



The same chain of arguments leads to critical exponents 
77 = 4 — d, 8 = d — 1, coincident with those of the d- 
dimensional Gross-Neveu model in the l&rge-Nf limit 
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[l9| . The nature of the transition predicted by SDEs 
thus appears sensitive to the order of the limits g 2 — > oo, 
Nf - N fc . 

We have been motivated by these considerations to 
attempt an EoS fit to the data of Fig. [4] of the form 

m = A[(N f - N fc ) + CL-^iWy + B{^) s , (11) 

where L is the linear extent of the system and the 
exponent v is given by the hyperscaling relation v = 
(S + l)/d(S—p). The term proportional to C accounts for 
finite volume corrections. Our best fit to data with Nf £ 
[4, 8] yields N fc = 6.89(2), S = 6.90(3), p = 4.23(2), with 
X 2 /dof=129/15. 

Besides the relatively poor quality, and significant dis- 
agreement with Eq. ([6]), the fit should not be taken too 
seriously for two further reasons, namely the combination 
of "exact" HMC data with HMD containing a systematic 
error due to At ^ 0, and the as yet unquantified errors 
due to the identification of gf im with g^ ax . Nonetheless 
the disparity of the value of the exponent S with the 
SDE prediction is striking. Moreover, if we plot the fit- 
Si i i i i i i i i ' 
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ted values for 6 and v together with those from EoS fits 
to data from Nf < Nf c compiled from Refs. [H, fl5| ■ 
as shown in Fig. [5j we see that the estimate for S(Nf c ) 
lies within a general trend of S increasing with Nf , with 
no evidence for non-commutativity of the limits g 2 — > oo 
and Nf — » Nf c (there is some indication for a jump from 
v 0.5 to v » 1 as Nf increases from 6 to Nf c , which 
needs to be confirmed in a more refined study). In any 
case, the fitted values of 6 lie well above those predicted 
in the SDE approach. 

In summary, for the first time we have been able using 
lattice Monte Carlo simulation to study a chiral phase 
transition as the number of fermion flavors Nf is var- 
ied continuously. We find a continuous phase transition 
in the strong coupling limit at a critical flavor number 
Nf c — 6.6(1), in qualitative, but not quantitative agree- 
ment with the predictions of an analysis using Schwinger- 
Dyson equations. In particular, the critical exponents are 
not those either of a conformal phase transition or the 3d 
Gross-Neveu model. It is plausible that our result may 
inform estimates of the corresponding critical number of 
flavors for chiral symmetry breaking in QED3, where di- 
rect lattice simulations are hampered by a large separa- 
tion of scales. Note that a recent perturbative analysis 
of RG flow equations in the large- Nf limit of QED 3 pre- 
dicts Nf c = 6 [20]. In future it will be interesting to 
check whether the same universal features of the strong 
coupling limit emerge using the alternative lattice formu- 
lations of the Thirring model reviewed in [l3| • 
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